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The edge of chaos is analyzed in a spatiaUy extended system, modeled by the regularized 
long-wave equation, prior to the transition to permanent spatiotemporal chaos. In the presence 
of coexisting attractors, a chaotic saddle is born at the basin boundary due to a smooth-fractal 
metamorphosis. As a control parameter is varied, the chaotic transient evolves to well-developed 
transient turbulence via a cascade of fractal-fractal metamorphoses. The edge state responsible for 
the edge of chaos and the genesis of turbulence is an unstable travelling wave in the laboratory 
frame, corresponding to a saddle point lying at the basin boundary in the Fourier space. 
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There is a growing interest in edge states at the 
laminar-turbulent boundary which can improve our un- 
derstanding of the transition from turbulent to laminar 
flows in fluids and plasmas, as well as the precursors of 
turbulence [TJ |5j . Dynamical systems description of the 
transition from temporally to spatiotemporally chaotic 
attractors, based on the regularized long-wave equation 
(RLWE) , provides a simple model to acquire in-depth in- 
sights of the laminar-turbulence transition . In this 
paper we use the RLWE to study the nonlinear dynamics 
of a spatially extended system prior to the onset of per- 
manent spatiotemporal chaos. The aims are threefold. 
First, we establish the link between the concept of edge 
of chaos (EOC) at the boundary of laminar-turbulent 
transition and the concept of chaotic saddle at the basin 
boundary of coexisting attractors. Second, we show that 
a chaotic saddle is born in a smooth-fractal metamorpho- 
sis which evolves to well- developed transient turbulence 
via fractal- fractal metamorphoses. Third, we elucidate 
the role of the edge state at the basin boundary of coex- 
isting attractors and at the boundary of pseudo basins of 
coexisting chaotic saddle and attractor before the onset 
of permanent spatiotemporal chaos, and at the boundary 
of pseudo basins of coexisting chaotic saddles/attractors 
after the onset of permanent spatiotemporal chaos. 

The driven-damped regularized long-wave equation is 
given by [31 U 

dtu + cdxU + fudxU + adtxxU = —uu—esin^Kx — flt) (1) 

where e is the driver amplitude, c = 1, / = —6, a — 
-0.28711, 1/ = 0.1, K = 1 and = 0.65. We impose 
periodic boundary conditions u{x,t) — u{x + 2n,t) and 
solve Eq. ([T]) numerically using a pseudospectral method 
by expanding the wave variable u{x,t) in a Fourier se- 
ries u{x,t) = J2k=-N''^k{t)exp{ikx). We set TV = 32 
[4j. Since u{x,t) is a real function, only fc > need to 
be considered, and after dealiasing 20 complex Fourier 



modes are kept, thus the phase space has dimension 40. 
In the absence of driving-dissipation (e = = 0) or when 
driving-dissipation is relatively weak (e,!^ < 1), Eq. ([T]) 
admits a steady wave solution in the form of a solitary 
traveling wave [S] . If we keep all parameters in Eq. ([T]) 
fixed and only vary e, the steady wave solution of Eq. 
([T]) eventually becomes unstable and undergoes a diver- 
sity of bifurcations, giving rise to a wealth of dynamical 
phenomena. 

At e = 0.199, just before the onset of permanent spa- 
tiotemporal chaos, the solutions of Eq. ([T]) exhibit the 
characteristics of edge of chaos. A technique to detect 
the edge of chaos is to compute the lifetime of initial con- 
ditions in some region of the phase space [2 [5] , defined 
as the time a trajectory takes to converge to the laminar 
attractor. Figure [TJ a) shows the lifetime landscape in a 
two-dimensional projection of the phase space. The red 
regions indicate short lifetime, and corresponds to initial 
conditions whose trajectories do not show the features 
of transient turbulence (governed by a spatiotemporally 
chaotic saddle STCS J^J) and converge quickly to the lam- 
inar attractor (spatially regular and temporally chaotic 
attractor). On the other hand, the light blue regions cor- 
respond to initial conditions whose temporal evolution 
displays long chaotic transients before converging to the 
laminar attractor. The stable manifold of STCS is well 
approximated by the regions of longer lifetime. The edge 
of chaos is the boundary dividing the two regions of life- 
time in Fig. [ija). 

The cross in Fig. [ija) marks the position of the edge 
state (ES), which lies on the edge of chaos. The edge 
state is found through the bisection method [T]. By 
integrating many different initial conditions it is seen 
that the trajectories associated with turbulence have 
high level of energy bursts; here, energy is defined by 
E — fg^[u{x,t)'^ — aUx{x,t)^]dx/4:'!T . In contrast, the 
trajectories that converge quickly to the laminar attrac- 



2 





0.02 0.06 0.10 0.14 0.18 
lu,l 




FIG. 1. Edge of chaos is the boundary that separates two 
regions in the 2D projection of the phase space in (a), at e = 
0.199, showing the transient lifetime for turbulent trajectories 
to converge to a laminar attractor. Edge state is indicated 
by a black cross whose stable manifold is the edge of chaos. 
The blue (red) regions indicates long (short) lifetimes that 
correspond to the initial conditions that do (do not) exhibit 
transient turbulence before converging to a laminar attractor, 
as illustrated by the time series of energy E in (b), obtained 
by the bisection method for two different initial conditions. 



tor have low level of energy fluctuations. Beginning with 
two initial conditions us and ul, with short and long life- 
times, respectively, we integrate the condition given by 
the middle point of the path that connects both condi- 
tions, um — (ms+wl)/2, until it converges to the laminar 
attractor. We set the energy level Eq = 0.2 as a thresh- 
old to decide to which region of the phase space um be- 
longs. If the maximum energy along the trajectory of 



Mm is lower than Eq, um lies in the laminar pseudo-basin 
and at the next step we set ml = um- Otherwise, um 
belongs to the turbulent pseudo-basin, hence at the next 
step ug = Um- Repeating this procedure, we find pairs of 
conditions at both sides of the edge of chaos, arbitrarily 
close to each other. Figure [l|b) shows an example of two 
initial conditions determined by the bisection method, 
with the distance between them ||us — ml|| < 10""'^^. The 
red curve is the trajectory of laminar condition ug , and 
the light blue curve is the trajectory of turbulent condi- 
tion Ml. As the inset in Fig. [ijb) shows, both trajectories 
remain close to each other initially with E remaining al- 
most constant until t ~ 50. That part of the solutions 
corresponds to trajectories passing near the stable mani- 
fold (EOC) of ES, which is a saddle point in the Fourier 
space, with constant energy for a given control parame- 
ter. For t ^ 50, the 2 trajectories separate quickly. The 
laminar trajectory converges immediately to the laminar 
attractor, while the turbulent trajectory traverses first 
the vicinity of a chaotic saddle before converging to the 
laminar attractor at t ^ 4000. By applying systemat- 
ically the bisection method it is possible to find a long 
trajectory close to ES. 

At e = 0.199, STCS, which is irregular in space 
and chaotic in time, computed by the stagger-and-step 
method |6j has a Lyapunov spectrum with 14 posi- 
tive Lyapunov exponents and a Kaplan- Yorke dimension 
~ 36; ES is a saddle point, with one positive Lyapunov 
eigenvalue and 39 negative Lyapunov eigenvalues, whose 
stable manifold separates the two regions of pseudo- 
basins in Fig. [ij^a); the laminar attractor is spatially 
regular and temporally chaotic with one positive Lya- 
punov exponent and a Kaplan- Yorke dimension of ^ 22 

a- 

In order to probe the origin of the edge state and the 
genesis of transient turbulence related to the aforemen- 
tioned EOC, we construct a detailed bifurcation diagram 
in Fig. [2] for ii^ as a function of e for attractors, chaotic 
saddles and ES. We adopt a Poincare map by plotting 
a point every time the trajectory obtained from Eq. Q 
crosses the plane |M2(i)| = 0.1 with d\u2{t)\/dt > 0. The 
chaotic saddle is computed by the sprinkler method [7]. 
For e =0 to 0.25 we have identified four different attrac- 
tors: Ai, A2, A3 and A4. In the interval e = to 0.079 
Ai is a stable fixed point (thick blue line) with a con- 
stant energy for a given e, which loses its stability and 
is converted to a period- 1 limit cycle (thin blue line) via 
a Hopf bifurcation (HB) at e ~ 0.079. We observe three 
small energy jumps in the Ai branch, one of them is 
visible in Fig. [2] at e ~ 0.154. These jumps represent 
transitions from one period-1 attractor to another. Al- 
though strictly speaking these are different attractors, 
we refer to them as Ai along the paper because they oc- 
cupy roughly the same area in the phase space and their 
bifurcations do not affect our main analysis. The last 
period-1 limit cycle vanishes at e ~ 0.1925. A2 appears 
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via a tangent bifurcation (TB) at e ~ 0.09 when two fixed 
points, one stable (thick magenta hne) and one unstable 
(ES, dashed black line), are created. This unstable fixed 
point corresponds to the edge state that plays a funda- 
mental role in the genesis of transient turbulence and 
EOC seen in Figs. [l]and[2] At e - 0.125, A2 suffers a 
Hopf bifurcation and becomes a limit cycle of period- 1 
(thin magenta line). At e ~ 0.1297, A2 is bifurcated into 
a quasiperiodic attractor which loses its stability and van- 
ishes at e ^ 0.13235. The coexistence of attractors Ai 
and A2 in the interval e ~ 0.09 to 0.13235 implies the 
existence of two basins of attraction. We will show that 
the bifurcation of the basin boundary is responsible for 
the genesis of transient turbulence. A period-2 attrac- 
tor ^3 (red line) appears via a saddle-node bifurcation 
(SNB) at e ~ 0.1774, which undergoes a number of dif- 
ferent bifurcations as e increases, involving a transition 
to quasiperiodicity, period-doubling cascade, and unsta- 
ble dimension variability to temporal chaos at e ~ 0.1925 
[5], and at e 0.2 it loses its stability via an interior cri- 
sis (IC) leading to the onset of a spatiotemporally chaotic 
attractor A4 (green) [1]. Chian et al. [J] demonstrated 
that for e % 0.21 A^ is composed of a spatiotemporally 
chaotic saddle (STCS, light blue) which preexists as the 
transient turbulence prior to IC and a temporally chaotic 
saddle (TCS, grey) evolved from ^3. TCS turns into a 
temporally chaotic attractor ^3 at e ^ 0.22105 due to 
a boundary crisis (BC). As e increases further, A^ turns 
into a period-1 limit cycle via an inverse period-doubling 
cascade and becomes a stable fixed point (thick red line) 
via a Hopf bifurcation at e ~ 0.2308. At e 0.235, the 
stable fixed point Ai disappears in a tangent bifurcation 
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FIG. 2. Bifurcation diagram of -E as a function of e for edge 
state (ES), attractors (A1-A4), spatiotemporally chaotic sad- 
dle (STCS) and temporally chaotic saddle (TCS), showing 
Hopf bifurcation (HB), tangent bifurcation (TB), saddle- node 
bifurcation (SNB), interior crisis (IC) and boundary crisis 
(BC). 



(TB), along with ES. A3 and A^ coexist for e = 0.22105 
to 0.2308. EOC can be found for 0.11 < e < 0.2 where 
there is coexistence of transient turbulence (STCS), edge 
state (ES) and spatially regular laminar attractors (Ai, 

Next we investigate the origin of the edge state and its 
role in the genesis and evolution of STCS responsible for 
the transient turbulence. As mentioned earlier, when A2 
appears as a stable fixed point, an unstable fixed point 
appears simultaneously via TB; Ai coexists with A2 for 
0.09 ^ e < 0.13235. The basins of attraction are sepa- 
rated by a boundary. ES is a saddle structure that lies at 
the basin boundary. We applied the bisection method [T] 
to detect ES that separates Ax and A2, and discovered 
that ES is the unstable fixed point born at the tangent 
bifurcation TB (e - 0.09), shown in Fig. [2] The edge 
state corresponds to an unstable travelling wave moving 
with the driver speed in the laboratory frame. 

Figure |3] shows the basins of attraction of Ai (blue) 
and A2 (magenta) for e = 0.095, 0.111, and 0.130, re- 
spectively, along with the edge state (cross). The sta- 
ble manifold of ES is the basin boundary separating 2 
coexisting attractors. For e < cm 0.11 , the basin 
boundary is smooth as shown by Fig. [sj^a) for e = 0.095. 
At e = Em, a chaotic saddle STCS associated with Smale 
horseshoe structures is born when the basin boundary be- 
comes fractal as the result of infinite homoclinic crossings 
between the stable and unstable manifolds of ES due to 
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FIG. 3. Basins of attraction for the coexisting attractors A\ 
(blue) and A2 (magenta) at e = 0.095 (a), 0.111 (b), 0.130 
(c). The black cross denotes the edge state. 
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a smooth- fractal metamorphosis [5], as seen in Fig. [sjb). 
As e increases for e > em the basin boundaries become 
increasingly complex due to a cascade of fractal-fractal 
metamorphoses, as shown in Fig. [Sjc). 

To measure the fractal dimension of the basin bound- 
ary for 0.09 ^ e ^ 0.13235, where Ai and A2 coexist, we 
compute the uncertainty exponent a = D — d, where D 
is the dimension of the phase space and d is the fractal 
dimension of the basin boundary, defined by the uncer- 
tainty fraction [TU] f{5) ^ (5", where 5 is the radius of the 
sphere of uncertainty. The computed a for this interval 
is shown by the magenta circles in Fig. |4ja) where we 
use a two-dimensional phase space projection |u2|) 
with D = 2. Figure |4|a) shows that for e ^ 0.11 a is 
w 1 implying a smooth basin boundary. For e }^ 0.11 
a becomes < 1 implying a fractal basin boundary. As e 
increases further, a steadily decreases implying the oc- 
currence of a cascade of fractal-fractal metamorphoses 
[9], in agreement with changes in the basin boundaries 
seen in Fig. [3] 

By noticing that the structure of the basins of attrac- 
tion of Ai and A2 just before the disappearance of A2 
is very similar to the lifetime function for Ai just after 
the disappearance of A2 in the same region of the phase 
space, we use this information to compute the fractal di- 
mension of the pseudo-basin boundary for 0.14 < e < 
0.19. We define the lifetime difference for 2 initial condi- 
tions separated by 6 as AT{r) — \T{r + 5) — T(r)|, where 
T(r) is the time an initial condition r takes to converge 
to Ai. We classify an initial condition r as uncertain if 
AT(r) > Ath, where Ath is the time-difference thresh- 
old. The uncertainty exponent for each Ath is given by 
the slope of the line fitted by a linear regression. The val- 
ues of f{6) calculated for the interval (ii) are indicated 
by the blue squares in Fig. |4]ja). 

To compute the fractal dimension of the pseudo-basin 
boundary for e — 0.199, we use the expression for the 
upper bound of the uncertainty exponent as a function 
of the mean lifetime of the chaotic saddle t [H |H] and 
the maximum Lyapunov exponent A,nax, a < l/(rA„iax) 
[llj . Chian et al. [4] accurately obtained Amax = 0.1405 
and T = 353.3 (in units of the driver period) for e — 
0.199, thus a calculated from the upper bound formula 
is 0.002018, denoted by the green triangle in Fig. |4]^a). 

To quantify the degree of spatial disorder of the chaotic 
saddle created at e ^ 0.11, we compute the time-average 
of the Fourier power spectral Shanon entropy [4] of 
STCS, given by SAit) = - X]fc=i Inpfe(i), where 

Pk{t) = lukit)]"^ / Ylk=i l^fc(*)P- Figure Qa) shows that 
the degree of spatial disorder increases with e until it 
reaches the maximum value near e ~ 0.199. 

The increase of the degree of spatial disorder with in- 
creasing driver amplitude is accompanied by an increase 
of temporal chaos. A way to estimate the number of 
positive Lyapunov exponents is to compute a large num- 
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FIG. 4. Variation with e of: a) the uncertainty exponent a 
(magenta circles, blue squares, green triangle) and the time- 
average power spectral entropy {SA{t)) (cross), b) the num- 
ber of positive Lyapunov exponents A'^^ (blue circles) and 
the maximum Lyapunov exponent Amax (red triangles). eM 
indicates the genesis of a chaotic saddle (STSC). 

ber of finite-time Lyapunov spectra of STCS and taking 
the average [12]. Figure |4][b) shows that the number of 
positive Lyapunov exponents (blue circles) increases 
steadily with increasing e, reaching its maximum value 
of iV+ = 14 at e ^ 0.199. Similar behavior is observed 
for the maximum Lyapunov exponent Amax (red trian- 
gle), shown in Fig. |4|^b). Figure [4] provides a consis- 
tent overview of the genesis and evolution of the tran- 
sient turbulence showing that the degree of complexity 
of transient turbulence (STCS) increases as e increases 
and evolves to a well-developed transient turbulence be- 
fore the transition to permanent spatiotemporal chaos. 

In conclusion, we have demonstrated that prior to the 
onset of permanent spatiotemporal chaos the regularized 
long- wave equation exhibits the behavior of edge of chaos 
(EOC), whereby a trajectory traverses a transient turbu- 
lent state before converging to a laminar state. The edge 
state responsible for the EOC and the genesis of turbu- 
lence was identified and a sequence of metamorphoses of 
the EOC was shown to be responsible for the appear- 
ance of a chaotic saddle and its subsequent evolution to 
a well-developed transient turbulence. Our results pro- 
vide a much clearer picture of the origin of turbulence 
in the regularized long-wave equation, which has been 
extensively studied as a general model of transition to 
spatiotemporal chaos [3H3 IH| ■ 
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